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We discuss under what conditions Brownian transport pro- 
cesses can display negative resistance. We prove it cannot 
occur on ID spaces like the circle or the line. We construct 
an entropic ratchet: an explicit two-dimensional model, and 
its collapse onto a branched ID backbone, showing negative 
resistance and rectification as a consequence of a geometric 
symmetry breaking. We derive an accurate numerical method 
for solving our 2D model. Finally, we discuss analogies and 
relevance to biological ion channels, in particular for channel 
inactivation and blocking. 



If a current arises as a result of some force, then it 
will flow "downhill" , in the direction in which it dis- 
sipates energy into heat; thus, the resistance is always 
positive. The diminishing of a current as the driving 
force becomes stronger is called negative incremental re- 
sistance, or just plain "negative resistance" (NR); devices 
that display NR exist and have important technological 
applications; these devices are, usually, also rectifiers. If 
energy is provided through suitable bias voltages, so as to 
displace the steady state of operation to the NR region, 
then the negative slope of the current can give rise to in- 
teresting instabilities like relaxation oscillations A 
typical such device from electronics is the tunnel diode. 
Much more importantly, the very existance of our ner- 
vous systems depends upon the ability to generate action 
potentials to transmit nerve impulses along axons; this 
phenomenon requires at least one NR rectifying device, 
which is known to be the sodium channel |j]-{7|. 

In this Letter, we will show how NR and rectification 
can arise in a Brownian transport process, through a 
mechanism we will call entropic trapping. Because of the 
purely geometric nature of this mechanism, we can call 
our model an entropic ratchet. 

We will study a Brownian particle in a periodic poten- 
tial, embedded in an equilibrium, constant temperature 
bath, and subject to a single force trying to advect it. 
This type of process is described by a Langevin equation 
of the form 

x = f (x) + £(t) 

where x belongs to some euclidean space, £ is Gaussian 
white noise satisfying (£i(t)£j(s)) — 2kT5ij5(t — s), and 
f is of the form -Vl^(x) + F with F constant, i.e., f is 
a vector field independent of time. We will assume V to 
be periodic along the direction of F. 

Why would it be interesting to search for NR in trans- 
port processes of this form? First of all, these processes 



have very wide applications |p|. Second, transport in 
symmetry broken potentials ( "ratchets" ) |8|-|l4| has been 
shown to be analogous to conduction through electronic 
diodes, capable of rectification; however, the mechanism 
through which a tunnel diode provides NR is intrinsi- 
cally quantum mechanical; hence, it would be interesting 
to provide classical analog. Finally, ratchet potentials in 
one dimension can be shown not to have NR, as we will 
now do. 

In one dimension, an equation of the form 

x = f(x) + £(t), (£(<)£(*)) - 2kT6(t - s) 
has an associated stationary Fokker-Planck equation 
d t P{x,t) + d x J(x,t) =0 

JO, t) = fP(x, t) - ktd x P(x, t) 

whose steady state, for periodic f(x), can always be 
solved in quadratures. If / has a zero spatial average, 
then it is the gradient of a periodic potential, and detailed 
balance and Boltzmann weights hold. If / does not have 
a zero average, then it can be written as / = — d x V + F 
with V periodic, and the stationary state can be solved in 
double quadratures The Fokker-Planck probability 
current J as a function of F is given by: 

m . Sf^> (!) 

Q(F) = j j^- v "^ F "-'^ 2 - F ^''-'U x 'dx (2) 

where e^. T means exp(x/kT). Q is positive; moreover, 
8fQ is also positive, and satisfies kTdpQ < 2nQ. This 
inequality implies dJ/dF > and hence Brownian trans- 
port in one-dimensional periodic potentials (under a 
steady force and in white noise) cannot show NR. (We 
will show later that this only holds if the underlying space 
is topologically trivial). 

Now we will construct an explicit example in two di- 
mensions. First we illustrate the notion of entropic bar- 
riers, through a very simple example. Let us consider the 
following potential 

V(x,y) = x 2 (cos(y) + l.l) 

We see that the potential is identically zero on the y axis, 
and bigger than zero everywhere else. Hence there are no 
true energy barriers impeding motion of a Brownian par- 
ticle along the periodic direction y. However, the shape 
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of the potential around the y axis is also important at 
nonzero temperature. For any given y, a slice of the po- 
tential along that value of y is a parabola; however, as 
motion progresses along y this parabola opens and closes 
periodically. In the absence of an external force, a Brow- 
nian particle will spend more time around y = n than on 
the bottleneck y — 0, every now and then jumping one 
period up or down in y, as if there actually were an energy 
barrier; these are called entropic barriers because, unlike 
a true activation energy, the timescales they induce do 
not follow Arrhenius-Kramers laws. In other words, as 
will see, the violation of these laws is responsible for the 
non-monotonic behavior of the current as a function of 
the applied force F. 

We will now break the parity symmetry along the y 
axis; if we loosely call any parity broken potential a 
"ratchet" , we can say the following potential is an en- 
tropic ratchet: 

V(x,y) =x 2 (cos(y + lncoshj;) + l.l)/2 (3) 

where In cosh a; is just an easy way to make a function 
that is both even and linear in x for large x. The equipo- 
tentials are now symmetry broken and look like a her- 
ringbone pattern, see Fig. 1. Our dynamics will be given 
by 

x= -W + Fy + £(t) (4) 

where y is the unit vector along y, and the noise corre- 
lators are as before. 

So, if we apply a force F along y, if the force is positive 
the particle will move forward without problem. How- 
ever, if the force is negative, the particle will move back- 
wards, but every now and then it will get into a spine 
of the herringbone. It will progress upon the spine for a 
while, deeper the stronger the force, and then will have 
to go against the force in order to climb back up and 
get again on the backbone. Thus it will be locked for a 
while, because the energy required to climb back up is 
a honest activation energy, and the time required to do 
so obeys an Arrhenius law. But the probability that a 
particle on the backbone will get into the spine does not 
depend exponentially on F, and so the particle spends a 
larger proportion of time blocked away in the spines as 
F becomes more negative. The overall timescale to get 
into and out the spine is not Arrhenius-like, yet near the 
bottleneck the potential energy of points in the spine and 
in the backbone is the same. If it was, the time the par- 
ticle spends stuck in the spine would cancel out with the 
time it spends in the backbone, and the current would 
be a monotonic function of F. 

We have evolved numerically Equation [| to obtain the 
average speed of the particle as a function of F and kT. 
The mean speed equals the Fokker-Planck current times 
the period (27r) of the potential. The result of our simu- 
lations is shown in Figure 2. 

In performing a numerical integration of Equation ^, 
one encounters several numerical problems. The simplest 



method for numerical integration is the Euler method, 

x(t + At) = x(t) + f(x)At + \/2kTAtn 

where the r\ are random gaussian numbers with unit vari- 
ance. While the Euler method is first order in At for an 
ordinary differential equation, it is only one-half-order for 
a Langevin equation, a property characteristic of diffu- 
sion. This implies that extremely small timesteps have to 
be used for accurate integration. Several methods have 
been proposed for increasing the order of the integration, 
including stochastic versions of Runge-Kutta . How- 
ever, the spines in our potential are a stiff problem, being 
long and skinny, which would limit the applicability of a 
Runge-Kutta scheme even in the absence of noise. The 
reason is that the timescale has to be small in comparison 
with the relaxation time on the fastest direction, in this 
case the skinny direction normal to the axis of the spine. 
But then the interesting time evolution is that along the 
axis of the spines, which becomes painfully slow. 

We have devised a method to cope with this problem. 
The fundamental problem in developing a Runge-Kutta 
scheme is that somehow one is assuming analyticity both 
of the vector field and of the solution; the last one is just 
not there. But one can safely assume analyticity of the 
vector field alone. If we expand the latter to first order, 

f (x + Ax) » f (x) + Ax • (Vf)(x) 

We use the fact that f is curl-free (and thus Vf is a sym- 
metric tensor) to diagonalize it and rotate to its eigenba- 
sis. In this base, the problem locally becomes a cross- 
product of independent Ornstein-Uhlenbeck processes. 
Thus, the question is: if we know we currently are at 
position x at time t, what is the probability that we will 
be at position x' at time t'l The answer to this question 
is the Fokker-Planck propagator (Green function), which 
is known analytically for the Ornstein-Uhlenbeck process 
H . It is a Gaussian, centered at the position x' that 
would be the solution for the deterministic case at time 
t' , and with a width which has been changed because of 
the compression/expansion due to /' = df / dx. Thus it is 
a trivial matter to generate a new value x' with the cor- 
rect probability distribution for an Ornstein-Uhlenbeck 
process; our numerical method then reads 

x' = {x- x c )e-J' At + x c + ((fcT//')(l - e- 2/ ' At )) ^ r, 

x c = x- f/f, At = t'-t 

for each eigendirection of the Hessian. The new coordi- 
nates are then rotated back to the original frame. The 
advantages of this method over stochastic Runge-Kutta 
are two. First, it solves exactly the Ornstein-Uhlenbeck 
process, by construction, and hence any linear problem. 
Second, it can handle stiff problems more easily; it will 
not lose accuracy if the timestep is larger than the re- 
laxation timescale of the fast direction, because it will 
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not overshoot and generate dynamical instabilities; one 
can concentrate on the more interesting slow timescale. 
The disadvantages are also two: analytic knowledge of 
the Hessian is required, and a matrix eigensystem calcu- 
lation has to be performed, together with two coordinate 
frame transformations. Thus our method rapidly loses 
ground to stochastic Runge-Kutta for high dimension- 
alities, unless the problem is quite stiff or the Hessian 
is sparse. For our stiff problem in two dimensions, this 
method is extremely well adapted; Fig 2 was generated 
in a few hours of cpu of a workstation and shows almost 
no trace of the residual noise typical of direct Langevin 
simulations. 

Back to our problem, we would like to get some more 
explicit understanding of the J(F) curve we have just 
obtained numerically. Let us consider a one-dimensional 
cartoon of the system: we can represent the backbone 
as a circle, with a single spine of length L attached to 
it, at an angle 9, as in Fig. 3a. There is no potential, 
just the steady force. At the vertex, we have to impose 
that the probability be continuous: the limit of P(x) as 
we approach the vertex from all three sides must be the 
same, P ver t- The current should be conserved at the 
vertex, so that probability does not accumulate. Then, 
the probability density on the circle is constant, and so is 
the current on the backbone: J = FP vert . The current 
on the spine should vanish, and hence the probability 
along the spine P(x) as exp(—Fx/ tan(0)). Therefore, 
the final current is 
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which has the right qualitative form (Fig 3b), except for 
a single detail. For F large and positive both cartoon and 
full case converge algebraically to unity. But for F large 
and negative, J — > exponentially, while in the full case 
it seems to die faster. This behavior can be understood 
if we recall that the spines are truly different: first, there 
is a potential along them; second, they are not of length 
L, but rather arbitrarily deep. The center of the spines 
lies approximately at y = ir — In cosh x , and the potential 
there is V = 0.05 * x 2 . The current in this case equals 



yjl + ^nh 2 (x)e-°- 05x2+Flncos ^dx 



(6) 



Thus, for a large and negative force, the particle entering 
the spine will find a stable fixed point at x as 10F. Thus 
the effective length of the spine increases as F becomes 
more negative; the activation barrier is the product of 
this effective length and F, and hence quadratic in F 
rather than linear. So we have shown that NR is actually 
possible in one dimension, but only if the topology of the 
space is more complicated than just a circle. 

We have presented an explicit example of a Brownian 
transport process showing NR, The NR region is gen- 
erated through an essentially entropic process. There 



are no obstacles, no energy barriers to the motion of the 
particle along the y axis. There is just a finite proba- 
bility of exploring space a bit out of the center, and end 
up blocked inside a spine. Since we are in a thermal, 
rather than quantum, situation, our system has much 
more similarity to biological ion channels than to elec- 
tronics. There are two well known instances in channels 
where NR is observed. The first one is in channels that 
can have an inactive state; this is the case of the Na 
channel, but some K channels also have this property; 
this is extremely important biologically, because the NR 
resulting from inactivation is essential to the regenera- 
tion of action potentials. The second one is the blocking 
of channels through large ions; this is important exper- 
imentally, because in order to assess properties of new 
channels, biophysicists will test for changes in behavior 
when the channel is "poisoned" with various compounds 
of known effect on known channels. The standard poison 
arsenal includes several large ions, that can get partially 
into the channel and block it; for example, tetraethylam- 
monium is used for Na channels jl6| and Mg jlj, Cs 
]l8[ or polyamines for K channels. It is worth not- 
ing that the I-V curves of such channels look extremely 
similar to that of our model, and, furthermore, several of 
the experimental measurements [p 161c| show the faster- 
than-exponential decay of our model, while standard the- 
oretical models with fixed barriers Jl7| , p0] | show exponen- 
tial decays; this discrepancy can only be solved through 
models having barriers that depend on the field. This 
might mean that the large ion buries itself deeper and 
deeper into the crevice of the channel, getting more and 
more stuck and having to climb a larger distance against 
the potential to get out, just as for our spines. There 
are also similarities with the case of channel inactiva- 
tion. These channels are hypothesized to have three 
states: closed, open, and inactive. The transition be- 
tween open and inactive states has been modeled with 
a "ball and thread" mechanism 0,^,^l|: some mobile, 
charged, globular piece of the channel loosely attached 
through a long thread can get stuck in the mouth of the 
channel and block it; pretty much like a large ion blocker. 
Our model does not include a closed state, since the cen- 
tral backbone is always "open"; thus we do not see the 
current going to zero as F becomes large and positive. 
The closed state can be bypassed by looking at the peak 
sodium current from a pulse rather than the stationary 
current p2j , and in this case the observed peak I-V curve 
is qualitatively similar to what we observe. 

We wish to thank Albert Libchaber, Sanderman Simon 
and Gustavo Stolovitzky for many discussions. 
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FIG. 1. The equipotentials of the model. 

FIG. 2. J{F) at kT = 1 as computed from Langevin 
simulations. Each run lasted for 10 5 units of time; ten runs 
were done for each value of F. 

FIG. 3. (a) The phase space for the cartoon, (b) J(F) as 
given by Eq. 4 
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